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ABSTRACT 

o' 

Context. The ESO-Spitzer Imaging extragalactic Survey (ESIS) is the optical follow up of the Spitzer Wide-area Infra- 
-^ ' Red Extragalactic survey (SWIRE) in the ELAIS-S1 region of the sky. 

Aims. In the era of observational cosmology, the main efforts are focused on the study of galaxy evolution and its 
environmental dependence. Wide area, multiwavelength, extragalactic surveys are needed in order to probe sufficiently 
large volumes, minimize cosmic variance and find significant numbers of rare objects. 

Methods. We present VIMOS I and z band imaging belonging to the ESIS survey. A total of ~ 4 deg 2 were targeted in 
/ and ~ 1 deg 2 in z. Accurate data processing includes removal of fringing, and mosaicking of the complex observing 
pattern. Completeness levels and photometric uncertainties are estimated through simulations. The multi-wavelength 
data available in the area are exploited to identify high-redshift galaxies, using the IR-peak technique. 
Results. More than 300000 galaxies have been detected in the / band and ~ 50000 in the z band. Object coordinates are 
defined within an uncertainty of ~ 0.2 arcsec r.m.s., with respect to GSC 2.2. We reach a 90% average completeness at 
23.1 and 22.5 mag (Vega) in the I and z bands, respectively. On the basis of IRAC colors, we identified galaxies having 
the 1.6/im stellar peak shifted to Z — 1 — 3. The new /, z band data provide reliable constraints to avoid low-redshift 
interlopers and reinforce this selection. Roughly 1000 galaxies between z = 2 — 3 were identified over the ESIS ~ 4 
deg 2 , at the SWIRE 5.8/im depth (25.8 /xjy at 3a). These are the best galaxy candidates to dominate the massive tail 
(M > 10 11 M Q ) of the 2 > 2 mass function. 

Key words. Surveys - Galaxies: evolution - Cosmology: observations - Galaxies: high-redshift - Galaxies: statistics 

1. Introduction mogeneous Universe detected right after the Big Bang into 

the clumpy, clustered structures we see at z = 0. 
In the recent years, the observational effort in cosmology 

has focused on the understanding of galaxy formation and With this respect, the critical parameters to be probed 

tracing the processes that transformed the smooth and ho- turned out to be baryonic mass and the epoch of its ma- 



jor assembly (e.g.lPozzettije t al. 20 Q7_ iRudnick et aT1 [2006, 



120031 : iDickinson etTaD l2003bt Brinchmann & Ellis 2000) 



Send offprint requests to: Stefano Berta, e-mail: In this respect, particularly controversial is the origin of 
berta@mpe.mpg.de, ste_atreb@yahoo.it the recently discovered population of massive (M > 10 1 1 



* Based on observations collected at the European Southern 
Observatory, Chile, ESO No. 168.A-0322(A). 



ESIS web page: http://www.astro.unipd.it/esis 



Mfp) galaxies at redshift z > 1 . 5 — 2.0 (e.g. iBerta et al.l 
200 7aFlFontana et all [2 006. 20 01 iFranceschini et all 120061 : 



Bundy et al. 2005 UDaddi et al.l 12004. ICimatti et ail l2002) 



Berta S., et al.: ESIS-VIMOS 



This puzzle has been approached with two opposite 
"philosophies" : very-deep, pencil-beam surveys, or wide- 
area, shallower surveys. 

The first approach has mostly relied on space-based 
and 8m telescopes observations, that provided the deepest 
and clearest views of the distant Uni verse so far. Some ex- 
amples, are the Hubble Deep Fields (Williams et al. 1996, 
2000), the GOODS survey ( Dickins onet alj)2003al) and the 
Hubble Ultra Deep Field (JBeckwith et all 12003). pushing 
the detection limit to the 30th magnitude. 

Due to their small areas, these surveys are affected 
by non negligible cosmic variance, which can be a sig- 
nificant source of uncertainty in the study of galaxy sta- 
tistical properties, eve n at high redshift (see for example 
ISomerville et~alll2004h . 

Moreover the search for rare objects, such as massive 
high-rcdshift galaxies, requires large volumes to be sam- 
pled, in orde r to detect t hem in representative numbers. 
For example, iBerta et al.l (|2007af ) find a number density of 
- 3 x 10" 5 [\4 Mpc -3 ] for M > 1.6 x 10 11 M galaxies at 
z = 2 — 3, which translates in only a few hundreds sources 
per square degree. 

In overcoming these effects, wide fiel d surveys, such as 
COSMOS (2 dee 2 IScoville et al.l [2007bl). VVDS (16 deg 2 
iLe Fevre et al.ll2004h. CFH TLS (4 to 410 deg 2 at different 
depths, e.g. iNuiiten et al.ll2005l ) play a crucial role. 

In this framework, the Spitzer Wide-area Infra-Red 
Extragalactic survey (SWIRE, 49 deg 2 ) is quite unique. The 
availability of Spitzer IRAC (3.6, 4.5, 5.8, 8.0 /j,m) data 
samples the restframe near-IR emission of galaxies up to 
redshift z = 3 at least, directly probing the stellar mass as- 
sembly of galaxies in the distant Universe. Despite SWIRE 
being a rather shallow survey (reaching 5er depths of 3.7 
and 43 fiJy at 3.6 and 5.8 fim, res pectively) , it detect s 
2 x 10 11 M galaxies up to z = 3 (jBerta et al.ll2007al) . 
MIPS (mainly 24/zm) data secure the identification of star 
forming systems in the same redshift range. 

An extensive multi-wavelength coverage is required. 
When dealing with large datasets over such a wide area, 
a spectroscopic characterization of the survey is beyond 
the multiplexing capabilities of current instrumentation, 
and photometric investigation is clearly the only viable ap- 
proach to be carried out in a reasonable amount of time. 

The ESO-Spitzer wide-area Imaging Survey (ESIS) is 
an ESO Large Programme (P.I. Alberto Franceschini) , tar- 
geting the SWIRE ELAIS-S1 field. This region includes the 
absolute minimum of the Galactic 100 /an emission in the 
Southern sky ( 0.37 [MJy/sr ]). 

In paper-I (jBerta et alj 120061 ) we presented the ESIS 
WFI/2.2m observations in the B, V, R bands, illustrating 
also the potential of panchromatic studies of galaxies. Here 
VIMOS / and z band observations are presented, complet- 
ing the optical ESIS survey. 

The paper is organized as follows. Section [2] presents the 
survey and the observing strategy. Section [3] and Appendix 
|A1 describe the main steps in the data reduction, including 
astrometric, photometric calibration and mosaicking. The 
quality of data products is tested in Sect. [H while the con- 
tents of the released catalogs are described in Sect. [5) In 
Sect. |5] we show how VIMOS / and z band data help in 
the selection and in constraining the spectral energy dis- 
tributions (SEDs) of high-redshift galaxies. Finally, Sect. [7] 
summarizes our findings. 



2. Observations 

The ESIS VIMOS project covers a total area of ~4 deg 2 
in the / band and ~1 deg 2 in the z band. A total of 116 
hours have been allocated, during periods P71-P73 (from 
April 2003 to September 2004); observations were carried 
over also in period P74 (from October 2004 to March 2005). 

Figure [1] shows the position of the surveyed area, with 
respect to other available observations. VIMOS fields are 
plotted as shaded areas (in blue/dark I band only and 
in red/light I + z pointi ngs). Multiwavelengt h imaging in- 
cludes Spitzer /SWIR E (lLonsdale et al.l 120031 l2004li. opti- 
cal E SIS WFI (IBerta et al.ll2006h. near-IR JK (iDias et all 
120071 ), GALEX (iMartin et al 1120051 : iBurgarella et al.ll2005lh 
XMM (jPuccetti et al.l I2006D. Chandra (Feruglio et al., 



sub.), ELAIS I SO (e.g. |R owan-Robin son et alj 12004 



ATCA 1.4 GHzjMiddclbcr 
Il999f ). See IBerta et alj ([ 200 



et al.l [20071 : iGruppioni et al.l 

1 -w^j for a summary and descrip- 
tion of all available data. The VIMOS z band observations 
are centered on the sub-area targeted by X-ray and JK 
observations. 



0.8 



RA (hours) 
0.7 ATC/ 0-6 



0.5 



<D 

u 

Ofl 
<U 

-a 

u 
Da 
Q 




10 

RA (degree) 

Fig. 1. The ESIS VIMOS survey in ELAIS-S1 and the mul- 
tiwavelength coverage of the area. The shaded area repre- 
sents the 49 ESIS VIMOS pointings in the / (blue/dark) 
and I + z (red/light) bands. O yer plotted are: t he ESIS 
WFI sur vey (~ 30' x 30' squares IBerta et al. ll2006D and the 
SWIRE ( Lonsdale et aD l200lll200l. XMM dPuccetti et al.1 



2006), ELAIS I SO (e.g. iRowan-R obinson et al. 2004), 



ATC A 1.4 GHz dMiddelberg et al.l 120071 : IGruppioni et al.l 
119991 ) areas. 



The VIMOS ([LeFevre et al.l l2003h array is made of 4 
CCDs, each one covering ~ 6.5' x 8' in the sky, separated 
by large gaps (~ 120 — 140 arcsec wide, depending of which 
CCDs are considered). The pixel scale is 0.205 arcsec/pixel. 

In order to cover a contiguous area and minimize over- 
heads, the observing strategy was set as follows: 
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1. the total area is divided into 37 I and 12 / 
(see Fig. p. 

2. each pointing is effectively made of 3 different Observing 
Blocks (OBs), named "A, B, C". Between pointings A 
and B one-gap- wide x,y offsets exist; the same holds 
between B and C pointings (see Fig. [2 left panel). 

3. each one of these 3 OBs is made of 4 dithered exposures, 
in order to get rid of bad pixels and gap signatures (see 
Fig.[H right panel). 

4. care has been taken in order to observe at least one of 
the 3 OBs per pointing during a photometric night. In 
this way, it is possible to compute photometric offsets 
and correct the whole field to photometric conditions. 

5. the total number of individual science exposures is 588 
and 144 in the / and z band respectively. 

6. the total exposure time per OB is 600s in the I band 
and 1200s in the z band, causing saturation at / ~ 16 
[Vega]. 

7. as a result of offsetting (3 OBs per pointing), the fi- 
nal coverage is not homogeneous, but there are regions 
with a 600s coverage, others with 1200s (the majority) 
and finally some with 1800s coverage. These values re- 
fer to the / band; as far as the z band is concerned, the 
effective exposure times are 1200s, 2400s, 3600s. 

Figure [3] shows the average VIMOS I and z pass- 
bands, defined as filters convolved with CCD and optics re- 
sponse. Actually the VIMOS transmission depends slightly 
on which of the four quadrants is considered. Here here we 
adopt an average transmission. Tables 0] and \5\ (available 
in electronic version) list the observation logs of the ESIS 
VIMOS survey. 



z pointings package, developed in Padova and adapted to the VIMOS 
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Fig. 3. The VIMOS I and z passbands average transmis- 
sivity. 



3. Data reduction 

The ESIS VIMOS data were reduced within the IRAfS 
environment, exploiting the capabilities of the WFPDRED 



De-biasing and flat-fielding were performed in the stan- 
dard manner; sky flat field frames taken during each night 
were used. VIMOS I and z band images are affected by 
strong fringing, that needs to be corrected before photo- 
metric calibration and mosaicking. Appendix [X] describes 
the main causes of fringing and the procedure adopted to 
correct ESIS frames. Figure H] shows one VIMOS I band 
frame before (left panel) and after (right panel) subtract- 
ing the fringe pattern (central panel). 

3.1. Astrometric calibration 

The ESIS survey targets a large area in the sky (~4 deg 2 
in the I band), therefore it is necessary to accurately map 
the astrometric distortion of VIMOS, in order to correctly 
build wide-field mosaics. 

Th e survey includes observations of the IStone et al.l 
(|l999f ) astrometric fields D and E, obtained both with the 
WFI and VIMOS instruments, and aimed at properly map- 
ping pixel into celestial coordinates. 

We have chosen a TNX mapping (see lBerta et alj '2006). 
which combines a linear projection of the sky sphere onto 
a tangent plane and a polynomial function for distortions. 
In order to obtain the best pixels-sky mapping, we have ex- 
ploited the VIMOS astrometric observations and extended 
the resulting catalog with the WFI observations. 

The best results were obtained using 5th order distor- 
tion polynomials. Before mosaicking, the astrometric so- 
lution was applied to individual science images, and each 
fram e was re-centered on the GSC 2.2 (STScI & OaTO, 
2001) catalog. Section I4TT1 reports on the accuracy of ESIS 
VIMOS astrometry. 

3.2. Photometric zeropoints 



1 The package IRAF is distributed by the National Optical 
Astronomy Observatory which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



As part of the ESIS VIMOS observations, lLandoltl (l992) 
Rul49 or TPHE spectrophotometric standard fields were 
targeted during photometric nights, in order to perform a 
dedicated photometric calibration of ESIS science frames. 
These observation have been defined with four different ex- 
posures per observing block, including the main standard 
stars on each individual CCD of the VIMOS array. 

Unfortunately, these fields are not rich in stars and the 
VIMOS chips are relatively small, therefore only few ob- 
jects fall on the CCD, and the estimate of the zeropoints is 
not optimal. Nevertheless, after the reduction of standard 
fields and correction of single CCD frames for differences in 
gain, the / band zeropoints were computed. The resulting 
zeropoint zpi = 26.96 ± 0.07 (in units of electrons) is con- 
sistent with the ESC0 value (26.99 on average), computed 
on richer fields. 

As far as the z band is co ncerned, despite the fact that 
ESIS included observations of lLando lt (1992) fields, no ref- 
erence standard catalog is available in these fields, there- 
fore we relied on a different procedure to determine the 
zeropoint. 

The ELAIS-S1 field benefits from an extensive multi- 
wavelength follow up from the X-rays to the radio (see 
iBerta et a l. 2006). We have therefore selected stars with 



2 http://www.eso.org/observing/dfo/quality/ 
VIMOS/qc/zeropoints.html 
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Fig. 2. Scheme of the ESIS VIMOS observing strategy. Left panel: each pointing is made of three different OBs, shifted 
in the sky in order to fill gaps between CCDs. Right panel: each OB consist of four dithered exposures, for optimal bad 
pixels masking. 
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(a) Image before de-fringing 



(b) Fringe pattern 



(c) Image after de-fringing 



Fig. 4. Example of de-fringing process in the I band. 



available BVRJK+IRAC photometry in the ESIS field and 
performed SED fitting to estimate the expected z band 
magnitudes and determine the VIMOS zeropoint: 



ZP 



^expected unmeasured ■ 



(i) 



Source extraction has b een performed with SExtractor 
(jBertin fc Arnoutall996f ) on individual science frames ob- 
served during photometric nights and reduced as described 
above. Each image was transformed into e~/s, by using the 
proper gain for the given CCD. 

Point-like objects were selected on the basis of the de- 
pendence of their half-flux radius on instrumental magni- 
tude (see Fig. bottom left panel, and Sect.E]). This pro- 
cedure has the advantage of well identifying the magnitude 
range over which point-like sources are well defined. 

The z band extracted catalog was matched to the other 
wavelengths source lists with a simple closest-neighbor algo- 



rithm, using a 1 arcsec matching radius (jBerta et al.l l2006) . 
The final catalog contains ^50-100 objects per CCD. Since 
we are using a sub-area of the entire ESIS VIMOS survey 
(i.e. the region in the sky actually covered by JK data), 
these photometric images were taken during two epochs 
only, namely December 2003 and Aug. -Sept. 2004. 

The adopted stellar library is based on iPicklesI (|1998f ) 
stellar spectrgQ extended to A > 2.5/im with a Rayleigh- 
Jeans law. 

The observed spectral energy distributions (SEDs) of 
stellar objects were fitted taking into account optical and 
near-IR photometry up to the IRAC 4.5/im band. Best fit is 
sought by means of \ 2 minimization, among all templates 
in the stellar library and varying their normalization. Two 
examples of SED fits are shown in Fig. \E\ The 5.8/im flux 



http://www.ifa.hawaii.edu/users/pickles/AJP/hilib.html 
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is also shown (when available), even if it was not included 
in the estimation of y 2 . The 1,2,3 a uncertainty on the 
expected z band magnitude, as derived from the best fit, is 
computed as xlest + 1, 4, 9. 

The zeropoint for a given CCD is given by the median 
of all zeropoints (one per star, see Fig. [5]). We used a 3 
s.i.q.ro clipping (corresponding to a 2 a clipping for a nor- 
mal distribution). 

As zeropoint differences between different chips are 
mainly due to gain differences, the resulting zeropoints for 
the 4 VIMOS quadrants are consistent with each other. 
The z band zeropoint obtained in this way is (in units of 
electrons): 



ZP(Z band) = 25.74 ± 0.11 



(2) 



where the uncertainty is given by 4 x s.i.q.r.. 

Similarly, we have taken advantage of the selected stellar 
catalogs to newly compute the zeropoint in the / band, di- 
rectly from ESIS data, obtaining a zeropoint of 26.93 ±0.13 
mag (in units of electrons) , consistent with our previous and 
ESO's estimates. No significant difference in zeropoints be- 
tween the two epochs is detected. 

3.3. Photometric calibration 

As already pointed out, the observations of ESIS VIMOS 
fields are spread across several nights, and are characterized 
by different sky conditions. 

The requirement that at least one OB per pointing was 
observed in photometric conditions facilitates the task of 
absolute calibration of the whole dataset; we call these 
photometric images "reference frames" . We have identi- 
fied point-like objects on each science image, in common 
with any reference frame and then computed the zero- 
point shift due to non-photometric conditions: Amag = 
mag an y - mag re f. 

Hence, each individual frame has been corrected to pho- 
tometric conditions, normalized to unitary airmass and ex- 
posure time, and converted to electrons: 



/, 



phot 



t v -tnOAiAM-Kx+Amag) 



te 



x CONAD, 



(3) 



where AM is the airmass of the given frame, K\ is the 
atmospheric extinction coefficient for Paranal (K\ — 0.027, 
K z = 0.02), and CONAD is the gain in e"/ADU. 

No color corrections are applied, and final catalogs are 
given in calibrated VIMOS magnitudes (i.e. not trans- 
formed to any standard photometric system) . The effective 
wavelengths for the / and z passbands are 8140 and 9050 A 
respectively. The flux zeropoints and AB shifts^ for these 
magnitudes are: 



S(I, Jy) = 2385 x io - 4 '""^ 7 - Ve ^) 
S{z, Jy) = 2207 x io - 4 "™^- v ^) 



mag{I, AB) 
mag(z, AB) 



mag(I, Vega) + 0.4763 
mag(z, Vega) + 0.5604 



(4) 
(5) 

(6) 

(7) 



as derived convolving the spectrum of Vega with VIMOS 
passbands (see Fig. [3]) . 



4 semi inter-quartile range 

5 mag AB = -2.5 x log (S„[erg s^cm" 2 !^- 1 ]) - 48.58 



3.4. Mosaicking 

The astrometrically- and photometrically-calibrated frames 
were finally combined together into the final ESIS mosaics. 

In the / band due to the large surveyed area and due 
to image size limitations, nine different mosaics were pro- 
duced, in a 3 x 3 pattern covering the whole ESIS-VIMOS 
field. For the z band it was possible to assemble one single 
final mosaic, ~ 1 deg 2 wide. 

During co-addition, all frames belonging to a given mo- 
saic were registered to the same astrometric map. Different 
mosaics have different astrometric maps, using the mosaic 
center as tangent point for re-projection. Bad pixels and 
"defects" (e.g. optical reflections or crossing of satellites, 
see the top panels in Fig. [9]) were masked and excluded 
from the final mosaics. 

Because of the complex observing strategy, the effective 
exposure time varies across the field: the majority of the 
area has a 75% depth, and smaller regions with 25, 50 and 
100% depths exist. We have therefore built weight maps 
to be associated with each scientific image during catalog 
extraction. We also include a coverage flag in the extracted 
catalogs, describing the actual depth at the position of the 
detected objects. 

The entire / band field is shown in Fig. O where we 
have scaled down and combined the 9 sub- areas, for display 
purposes only. The z band mosaic is in Fig. [JJ Figure [8] 
includes exposure maps. Unfortunately gaps in RA between 
different pointings are left. This problem is mainly caused 
by shadowing of the E-W edges of VIMOS CCDs in ESIS 
images (see the bottom left panel in Fig. [9|). The width 
of the shade varies from frame to frame, and consequently 
so does the width of the residual gap in the final mosaics. 
Note that this is different from shadowing by the guiding 
star probe (bottom right panel). The actual reason for this 
loss of pixels is not clear. 



4. Quality of data products 

Observations were carried out under very different atmo- 
spheric conditions, over a ~ 2 years period of time. The 
gaussian FWHM of point-like sources is therefore not con- 
stant across the ESIS field, but there are regions with seeing 
< 0.8 arcsec and others with FWHM ;> 1.0 arcsec. The final 
catalogs were extracted on the mosaics without smoothing 
to the worst seeing in the field, in order to take advantage 
of good imaging quality, where available. 

We have tested how many sources would be lost by de- 
grading the images with good seeing to a 1 arcsec FWHM 
PSF, by convolving all individual frames with the appropri- 
ate kernel. As a result, roughly 4-5% of the sources is lost 
when smoothing the images, mainly at the faintest magni- 
tudes, after the turnover in the number counts (see Sect. 

El). 

Figure[lO]shows the trend of seeing as a function of mag- 
nitude (left panels). The large scatter reflects what men- 
tioned above. The deviation at bright magnitudes is due 
to incoming saturation effects. The saturation thresholds 
for the ESIS VIMOS survey turn out to be J ~ 16.5 and 
z ~ 16 mag [Vega], consistent with the expectations of 
ESO's exposure time calculatoiQ (J ~ 16.2 and z ~ 15.8 
for a point-like source with 1" seeing and A.M. = 1.3). 



http://www.eso.org/observing/etc/ 
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The right hand panels of Fig. [10] illustrate the distri- 
bution of seeing across the field. Pointlike sources are plot- 
ted in different grey levels (darker indicating better see- 
ing) and dot dimensions (smaller for better FWHM) . Three 
bins are considered: FWHM < 0.8, 0.8 < FWHM < 1.0, 
1.0 < FWHM < 1.2. 



4.1. Astrometric accuracy 

The coordinates of the objects detected on the final mo- 
saics have been compared to CSC 2.2 (STScI & OaTO, 
2001) sources, after re-centering the GSC catalog on the 
ESIS images. Figure [Tl] shows the results for one of the / 
band mosaics (left panels) and the z band final image (right 
panels). The r.m.s. of the distribution of coordinate differ- 
ences for the / and z bands is 0.195 and 0.199 arcsec in 
A(RA ■ cos[Dec\) and 0.176 arcsec in A(Dec). 

In the central panels in Fig.[TT]the coordinate differences 
between the ESIS and GSC catalogs are plotted against RA 
and Dec: no systematic trends are detected. 

As far as the / band is concerned, we have checked also 
the coordinate match between the 9 different mosaics, an- 
alyzing the differences in (RA,Dec) for GSC sources in the 



overlap regions (left bottom panel in Fig. [TTj) . The r.m.s. 
of the distribution of coordinate differences are 0.105 and 
0.087 arcsec in A(RA ■ cos[Dec]) and A(Dec) respectively. 

4.2. Photometric accuracy 

The photometric accuracy of the data has been tested with 
simulations, adding synthetic sources to the ESIS final mo- 
saics (with IRAF), adopting a 1" seeing, and performing 
source extraction as for science data. 

In each band, we produced several simulated images, 
by adding point-like, De Vaucouleurs or exponential disk 
objects at random positions on the science frames. In each 
case, a different image was produced per 0.25 mag bin, in 
the range 18-27 mag. A population of ~ 8000 sources/deg 2 
was added each time. 

Figure [12] compares the input magnitudes of the sim- 
ulated point-like sources and the values measured by 
SExtractor. For each of the two bands, the four panels re- 
fer to regions in the mosaics having effective depths < 50%, 
50 - 75% and > 75% (see Sect. [§). The solid lines repre- 
sent the median differences, the short-dashed lines are the 
1-er standard deviation and the long-dashed lines trace the 
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Fig. 6. The ESIS VIMOS final 7 band mosaic. 



semi-inter-quartile ranges (s.i.q.r.). Typical s.i.q.r. uncer- 
tainties at mag = 21.5, 23.5 (Vega) are ~ 0.05, 0.20 in the 
I band and ~ 0.08, 0.22 in the z band, on average on the 
whole area. 

Slightly larger uncertainties are found for De 
Vaucouleurs profiles and exponential disks: the corre- 
sponding s.i.q.r. are ~ 0.07, 0.22 at mag = 21.5, 23.5 
(Vega) in the I band and ~ 0.10, 0.27 in z. The usual 
systematic offset for the recovered magnitudes of De 
Vaucouleurs objects is present. Its value is ~ 0.2 — 0.3 
mag nitudes. This offset is well described in the literature 
(e.g. iFasano et al J H9980 : it is due to missing flux in the 
outskirts of De Vaucouleurs profiles and does not depend 
on the adopted source extraction code. 



4.3. Completeness 

By comparing the numbers of input and detected sources 
in simulations, we derive the detection rate in the ESIS 
VIMOS survey. 

Figure Q2] shows the results for point-like sources in the 
I (left) and z (right) bands, split by coverage again. The 
vertical dotted lines set the 90% and 95% completeness 
levels. 



On average, the 90% level is reached at magnitude 23.1 
and 22.5 (Vega) in the I and z band respectively. These 
values increase to 23.6 and 23.0, if only the regions with 
coverage > 75% are taken into account. 

As far as exponential disks and De Vaucouleurs profiles 
are concerned, the 90% detection rate is shifted to ~ 0.3 
and ~ 0.5 brighter magnitudes levels, respectively. 

5. Catalogs 

The total exact area covered by the I band mosaics is 4.11 
deg 2 . Excluding the gaps and areas around very bright 
stars, the catalog covers 3.89 deg 2 . We define the cover- 
age flag as the fraction / of nominal exposure time (900s 
for the I band, 1800s in z) spent on source; 25%, 64% and 
96% of the sources in the I band catalog lie in regions with 
/ > 0.75, 0.50, 0.25 respectively. 

The z band survey includes a 1.08 deg 2 area. After re- 
moval of regions not covered and big stars, the catalog in- 
cludes 0.98 deg 2 , 20%, 55% and 96% of which are above 
/ > 0.75, 0.50, 0.25. 

Catalogs we re extracted using SExtractor 

(|Bertin fc Arnoutall996f ) from the final mosaics, weighting 
the data with coverage maps. A total of 312929 and 57926 
sources were measured in the I and z band respectively. 
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Fig. 7. The ESIS VIMOS final z band mosaic. 
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Fig. 8. The ESIS VIMOS final / (left) and z band (right) exposure maps. Brighter regions correspond to deeper coverage. 
North is up, east is left; the two images have roughly the same scale. 
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Fig. 9. Examples of defects to be masked before mosaicking. Top panels show emission features: optical reflections and 
satellite trails. Bottom panels exemplify "vignetting" features: the unidentified shade causing residual gaps in the final 
mosaics (left) and the shadow of the guide star probe (right). For clarity sake we have inverted colors of the shades. 
North is up, east is left. 



Column Name 



Description 



Units 



1 


ID 


2 


RA 


3 


DEC 


4 


MAG.BEST 


5 


MAGERR_BEST 


6 


MAG.AUTO 


7 


MAGERR_AUTO 


8 


MAGJSOCOR 


9 


MAGERRJSOCOR 


in 


MAG.APER1 


n 


MAG.APER2 


12 


MAG.APER3 


13 


MAG_APER4 


14 


MAG.APER5 


15 


MAGERR_APER1 


16 


MAGERR.APER2 


17 


MAGERR.APER3 


18 


MAGERR.APER4 


19 


MAGERR.APER5 


20 


SN 


21 


KRON_RADIUS 


22 


FWHMJMAGE 


23 


CLASS_STAR 


24 


FLUX_RADIUS 


25 


FLAGS 


26 


IMAFLAGSJSO 



ESIS ID 

Right Ascension (J2000) 

Decimation (J2000) 

Best of MAG.AUTO and MAGJSOCOR 

RMS error for MAG_BEST 

Kron-like elliptical aperture magnitude 

RMS error for Kron-like magnitude 

Corrected isophotal magnitude 

RMS error for corrected isophotal magnitude 

Aperture magnitude (1.230 arcsec diam.) 

Aperture magnitude (2.460 arcsec diam.) 

Aperture magnitude (3.280 arcsec diam.) 

Aperture magnitude (4.715 arcsec diam.) 

Aperture magnitude (6.675 arcsec diam.) 

RMS error for APER1 magnitude 

RMS error for APER2 magnitude 

RMS error for APER3 magnitude 

RMS error for APER4 magnitude 

RMS error for APER5 magnitude 

Local S/N ratio 

Kron aperture in units of A or B 

FWHM assuming a gaussian core 

Star/Galaxy classifier output 

Half-flux radius 

SExtractor extraction flags 

Coverage flag 



[deg] 

[deg] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 



[pixel] 
[pixel] 



Table 1. Contents of ESIS VIMOS / and z band catalogs. 



The list and description of columns included in the re- 
leased catalogs are provided in Tab. [TJ All magnitudes are 



given in Vega units. Data are provided as electronic tables 
associated to this work. 
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Fig. 10. Actual FWHM of pointlike sources in the ESIS VIMOS I band (top panels) and z band (bottom panels) survey. 
The dependence on magnitude and the variation of FWHM across the field are shown. Grey levels and dot dimensions 
illustrate seeing properties, darker and smaller dots referring to better FWHM. Three bins are considered: FWHM < 0.8, 
0.8 < FWHM < 1.0, 1.0 < FWHM < 1.2. 



Figures [TJ] shows the ESIS VIMOS number counts, 
as split by effective depth (top panels) and into point- 
like and extended sources (bottom panels ). The latter 
are compared to data from t he COSMOS (Scoville et al.l 
2007al iTaniguchi etHI l2007f). Cx OODS-N (|Capak et all 



200C WPS (|McCracken et al.l [2 003). HD FN, H DFS, 
Willia m Herschel (WH DF |Metcalfe_et_alJ l200l and 
SDSS (|Yasuda et al.ll200"LT) surveys. The ESIS data are fairly 
consistent with the literature; discrepancies can be due to 
cosmic variance at the bright magnitudes, to a different 
criterion in galaxy/point-like sources separation at inter- 
mediate fluxes, and to incoming ESIS incompleteness at 
the faint end. We warn also that saturation effects might 
play a role at the brightest fluxes (see Sect. BJ. 

The data collected from the literature used the 
SExtractor's stellarity flag to distinguish these two classes, 
while here the ha lf-flux radius is exploited (see for example 
iBerta et al.l |2006|) . Figure [T5] shows the different selection 
of point-like objects for the z band, based on a simple cut 
in the stellarity flag (> 0.95 in this case, right panel) or 
based on the dependence of the half-flux radius on magni- 
tude (left panel) . Note again the effects of pixel saturation: 
the bright point-like objects deviate from the locus defined 
by half- flux radius. The SExtractor stellarity flag misses 



the brightest saturated point-like sources and tends to be 
contaminated by galaxies at intermediate fluxes (see also 
IBerta et al.ll2006l ). 

Table[2]reports the observed number counts of the ESIS 
VIMOS survey. 

6. High-redshift galaxies 

Large area, multiwavelength extragalactic surveys are de- 
signed to study the evolution of galaxies across cosmic time. 
SWIRE-ESIS covers the whole electromagnetic spectrum 
from the X-rays to radio wavelengths and is well suited to 
identify large samples of any class of objects reachable at 
its depth. 

Of particular interest is the recently discovered class of 
massive galaxies (M > 10 11 M©) above redshift z > 1.5, 
whose formation and evolution appears to be faster than 
what predicted by pure hierarchical s cenarios and appar 



ently showing downsizing effects (e.g. Cimatti et al. 



Daddi et all 12 004: 



Fontana et aTll200a iBundv et al. 
2006H 



2002; 



2005: 



Franccschini et al. 

Thanks to IRAC observations, SWIRE-ESIS directly 
probes the assembled stellar mass of galaxies up to red- 
shift z ~ 3. Due to its moderate depth, SWIRE can detect 
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Fig. 11. Astrometric accuracy of the ESIS VIMOS survey. Left panels: I ban d. Th e top and central plots show the 
displacement of ESIS coordinates with respect to the GSC 2.2 (STScI & OaTO. l200lh catalog. The bottom figure shows 
the relative astrometric accuracy between the 9 different / band sub-areas. Right panels: same plots for the z band. Only 
absolute astrometry is shown, because the z band covers 1 deg 2 and consists of one single final mosaic. Horizontal solid 
lines set the ±lcr levels of A(coords) . 
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Fig. 12. Photometric accuracy, derived by extracting synthetic point-like sources from ESIS VIMOS / band (left) and z 
band (right) mosaics. Different panels refer to different effective depths. Solid, short-dashed and long-dashed lines trace 
the median, the 1-cr and the s.i.q.r. levels respectively. 
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Fig. 13. Detection rate as a function of magnitude, from simultaed point-like sources. The vertical dotted lines set the 
90% and 95% completeness levels. 



only the most massive tail of the high-redshift g alaxy stellar 
mass f unction (e.g. M > 2 x 10 11 M Q at z ~ 3. iBerta et al.l 
l2007al) . 

Here we describe how the availability of / and z band 
data, when combined with the other multiwavelength obser- 
vations, provides good constraints in the selection of high-z 
massive galaxies. 

For the curren t purpose, we h a ve ma t ched t he SWIRE 
Spitzer catalog fLonsd aIeet~aLl 120031 I2004D to ESIS 
B, V, R, I and z (Berta et al. 2006, and this work), near- 
IR J, K s (iDias et al.l 12007ft . and GALEX GR43 data, us- 
ing a closest neig hbor criterion with a 1 arcsec radius (see 
iBerta et a l. 2006). We use total magnitudes in all bands. 
The SWIRE 3cr depths are 2.2, 3.2, 25.8, 22.6, 200 /xJy 
at 3.6, 4.5, 5.8, 8.0 and 24 /mi; the ESIS BVR 95% com- 
pleteness levels are reached at magnitudes B V ~ 25 and 
B ~ 24.5 (Vega); the J, K s catalogs are 95% complete at 
19.82, 18.73 mag (Vega) respectively. 

The ESIS VIMOS / band area includes 219877 SWIRE 
sources, detected in at least one IRAC band. Among these, 



http://galex.stsci.edu/GR4/ 



142776 (65%) have at least one optical detection (BVRIz), 
137206 are detected in the / band, 52620 in BVR (one 
band at least, over 1.5 deg 2 ), 35670 in z (over ~ 1 deg 2 ), 
and 30355 are detected in the GALEX DR4 ELAIS-S1 sur- 
vey. The z band area, covered by all bands (BVRIz, JK S , 
Spitzer, GALEX, X-rays and radio), includes 63770 SWIRE 
sources, 45873 (~ 72%) of which detected in at least one 
optical band. Table [3] summarizes these statistics. 



6.1. Selection of IR-peakers 

Several selection criteria have been defined to identify high 
redshift galaxies, exploiting optical and n ear-IR color s (e.g. 
Ly-break galaxies, LB Gs, S teidel et al. Il996t [[999; BX- 
BM, Adelbe rger e t al. 12004 distant red galaxies, DRGs, 
Franx et al . 120031 ; extremely red objects, EROs, Cimatti 
et al . 120021 ; distant obs cured galaxies, DOGs, Dey et al. 
120081 BzK, Daddi et al. l200l . LBGs, BX-BM's and DOGs 
classes consist of very faint objects that require [/-band 
and/or deeper optical observations than those available in 
the ELAIS-S1 field. EROs and DRGs in ELAIS-S1 are being 
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Fig. 14. Observed source number counts for the ESIS VIMOS survey, split into effective dep th bins (top panels 
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discussed in a forthcoming paper by Dias et al. (2008, in of high-z (massive) galaxies by identifying the rcstframe 

prep.). 1.6/xm stellar peak (Sawicki 2002; Simps on fc Eisenhardtl 

Recently, iBerta et all (|2007al) and lLonsdale et~aTl I1999T ) in the SEDs of galaxies, redshifted to z = 1.5 - 3.0. 
(2006) have shown that IRAC data allow the selection 
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Fig. 16. Color-color plots adopted in the selection of IR-peakers. Filled triangles and circles represent 5.8/im-peakers 
and 4.5/im-peakers respectively. Open circles are 4.5/im-peakers identified as low-redshift interlopers on the basis of 
(I — 3.6) ab- Small dots trace the general SWIRE-ESIS population; orange dots are general sources with no 8.0/mi 
detection (upper limits are used). Redshift tracks fo r different templates are plotted: M82 (dusty starburst) and M51 
(spiral) belong to the GraSil library (jSilva et al.lll998t ); the Syl and Sy2 templates represent Mrk231 and IRAS19254-7245 
(|Fritz et al.ll2006l : iBerta et al.ll2003l ). The black thick solid lines in the BzK diagram (bottom- right panel) represent three 
GraSil spiral models with no extinction (ages 0.2, 1.0, 2.0 Gyr) and one GraSil elliptical (age 3 Gyr). Redshift values are 
reported along tracks. The diagonal and vertical dashed lines in the left-hand diagrams represent the color conditions to 
identify low-redshift interlopers (see text). 



We called these objects IR-peakers (informally known also 
as "IR-bump g a laxies" ) . 

IBerta et all (|2007al) have presented the selection of IR- 
peakers and studied their stellar content. Here we refine and 
reinforce the selection using / band data. This allows us not 
only to avoid low-redshift interlopers, but also to extend the 
search for these sources over the whole ESIS-VIMOS area 
(~ 4 deg 2 ), instead of restricting it to the central square 
degree where optical BVR imaging is available to date. 

The main IR-peak selection is based in the detection of 
the 1.6/im peak in the IRAC domain. Two different classes 
of galaxies are identified, those for which the SED peaks 
in the 4.5/im (z = 1 — 2) and those in the 5.8/im band 
(2 = 2-3): 



S„(3.6) < 5„(4.5) > Si,(5.8) 

S„(3.6) < 5„(4.5) < S„(5.8) > S„(8.0) 



(8) 
(9) 



5.8/im-peakers include also sources not detected at 
8.0/im, but whose 8.0/im upper limit is still consistent with 



Eq. [9] This selection is shown in the top pa nels of Fig. 
[TBI in the classic IRA C color-color diagrams (|Lacv et al.l 
2004:; IStern et al.|[2005T ). Triangles and circles depict 5.8/im- 
and 4.5/im-peakers. Open symbols represent low-redshift 
interlopers, identified on the basis of optical-IR colors, as 
described below. For clarity's sake, we do not plot upper 
limit arrows for IR-peakers not detected at 8.0/im. Redshift 
tracks for a variety of templates are shown, up to z = 3 
(see labels on plot). Such a selection has also the advan- 
tage to avoid any contamination from AGNs, which are 
characterized by a restframe near-IR power-law-like (mono- 
tonic) SED. At the SWIRE 5.8/mi 3cx depth (25.8 /xJy), 973 
5.8/im-peak galaxies and 938 4.5/im-peakers were identified 
(see Tab. E]). 

As shown bv lBerta et al.l (|2007al ). unfortunately 4.5/im- 
peakers suffer a significant contamination from low redshift 
interlopers, when the 4.5/im flux is boosted by a strong 
3.3/tm PAH emission. On the other hand, in the case of 
5.8/im-peakers, two photometric bands sample the SED on 
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/ band 


z band 


mag. 


log TV 


logiV 


[Vega] 


Tot. Ext. 


Tot. Ext. 



14.0 


1.915 


0.711 


1.991 


0.610 


14.5 


2.164 


1.267 


2.343 


1.087 


15.0 


2.390 


1.410 


2.388 


1.724 


15.5 


2.497 


1.736 


2.565 


2.254 


16.0 


2.526 


1.736 


2.579 


2.223 


16.5 


2.655 


2.069 


2.656 


2.366 


17.0 


2.738 


2.216 


2.792 


2.597 


17.5 


2.870 


2.430 


2.964 


2.750 


18.0 


3.024 


2.675 


3.113 


2.937 


18.5 


3.181 


2.945 


3.261 


3.143 


19.0 


3.361 


3.188 


3.437 


3.357 


19.5 


3.513 


3.409 


3.580 


3.536 


20.0 


3.659 


3.623 


3.709 


3.687 


20.5 


3.795 


3.788 


3.891 


3.888 


21.0 


3.938 


3.938 


4.017 


4.017 


21.5 


4.087 


4.087 


4.135 


4.135 


22.0 


4.216 


4.216 


4.232 


4.232 


22.5 


4.338 


4.338 


4.264 


4.264 


23.0 


4.410 


4.410 


4.180 


4.180 


23.5 


4.409 


4.409 


3.951 


3.951 


24.0 


4.270 


4.270 


3.433 


3.433 


24.5 


3.901 


3.901 


- 


- 


25.0 


3.008 


3.008 


- 


- 


Table 2. ESIS VIMOS 


/ and z 


band number 


units of [deg -2 i 


Tiag 1 ]. 


No completeness corre 


applied. 











Class 



area 

_2i 



number 



[deg 



SWIRE (/-band area) 


4.11 


219877 


SW1RE+ BVRIz (one band) 


3.89 


142776 


SWIRE+ BVR (one band) 


1.5 


52620 


SWIRE+/ 


3.89 


137206 


SWIRE+z 


0.98 


35670 


SWIRE+GALEX 


4.11 


30355 


24/xm 


4.11 


12770 


SWIRE (z-band area) 


0.98 


63770 


SWIRE+ BVRIz (one band) 


0.98 


45873 


SWIRE+ BVR (one band) 


0.98 


39883 


SWIRE+/ 


0.98 


41337 


SWIRE+z 


0.98 


35670 


SWIRE+GALEX 


0.98 


8884 


24fim 


0.98 


3733 



4.5/im-peak (total) 4~Tl 938 

4.5/im-peak (interlopers) 4.11 353 

5.8/im-peak AAA 973 

Table 3. Multi-wavelength detection statistics in the 
ELAIS-S1 area. 



the blue side of the 1.6/im peak, therefore the conditions 
S„(3.6) < S V (A.5) < S„(5.8) provide a robust criterion to 
avoid interlopers. 

Optical data provide a significant help in breaking the 
degeneracy and identifying interlopers in the 4.5/im-peak 
sample. The bottom-left panel in Fig.HBlshows the position 
of IR-peakers in the (I — 3.6)^5 vs. (3.6 — 5.8)ab space, 
over the whole ESIS-VIMOS area. According to template 
redshift tracks, those sources having 



are likely to be z < 1 interlopers . This r esult is consistent 
with what found by iBerta et al.l l|2007al) . based on (K s — 
3.6) and (Re — 3.6) colors. 

Thanks to / band data, interlopers are then identified 
on IRAC color-color diagrams (open circles). Here, the di- 
agonal dashed line depicts the condition (4.5 — 8.0)ab = 
1.5 x ( 3.6 — 5.8)ab + 0.7, already found in IBerta et aD 
(|2007af ). and defined thanks to K s band data. A total of 
353 4.5/im-peakcrs at redshift z > 1 are left. 

The availability of z ban d data over ~ 1 de g 2 allows us 
to build the BzK diagram (|Daddi et alj 12004 ) and verify 
where IR-peakers lie in this color-space (bottom-right di- 
agram of Fig. [T6|) . Only those IR-peakers detected in the 
BzK+IRAG bands are taken into account, including 38 
5.8/im-peakers and 30 4.5/xm-peakers (excluding interlop- 
ers) . The BzK diagram confirms that the IR-peak selection 
effectively identifies objects above z ~ 1.4, with our sources 
lying in the star-forming locus. The majority of (7—3.6) < 2 
interlopers fall in the z < 1.4 region, with a few sources con- 
taminating the high-z area. Roughly 30% and 40% of BzK 
4.5/im- and 5.8/im-peak galaxies are detected by MIPS at 
24/xm at the SWIRE depth. 

None of the BzK -detected IR-peakers falls in the pas- 
sive z > 1.4 galaxy class. The median K s magnitude of 
the general BzK population detected in the ESIS area is 
20.14 [AB]. BzK galaxies with (z - K s ) > 2.5 have a me- 
dian K s — 20.65 [AB]. A passive BzK galaxy (with me- 
dian K s magnitude), characterized by (z — K s )ab > 2.5 
and (B — z)ab J> 2.7, would be detected at z > 23.1 and 
B > 25.8 [AB]. Unfortunately, this goes beyond the limit 
of the ESIS survey. 

Finally, in addition to the confirmation coming from the 
test based of the BzK criterion, it is worth to recall that 
independent s tudies of IR-peaker s based on optical Keck 
spectrosc opy (IBerta et al. 2007bl) and IRS mid-IR spec- 
troscopy (Wccd man et al.l 120061 ) showed that the IRAC- 



(I - 3.6) < 2 [AB] 



(10) 



based broad-band selection successfully identifies galaxies 
above z ~ 1.5. 



6.2. Photometric redshift 

When dealing with a very large number of sources, as in 
the case of our wide area survey, the estimate of galaxy 
redshifts is mainly based on photometric techniques. 

Let us check the advantage of introducing VIMOS / and 
z bands in the estimate of photometric redshifts for the 
general galaxy population in the SWIRE-ESIS survey, and 
explore whether the data in hand are sufficient to produce 
reliable photometric redshifts in those areas not covered by 
other optical (e.g. BVR) observations. 

To this aim, we have exploited the capabilities of 
the newly public photometric redshift code EAZY (G. 
Brammer, P. G. van Dokkum & P. Coppi, Ap J submitted), 
which combines a maximum likelihood (ML) algorithm and 
Bayesian constraints on the expected redshift distribution 
of galaxies. Bayesian techniques in photometric redshift 
estimates have been thoroughly studied and ap plied to 
large datasets in the past with good results (e.g. iBenftezI 
2000). An interesting comparison bet ween Bayesian and 
pure M L techniques was performed by Hil debrandt et all 
(2008), who show the advantage of the former, especially 
in avoiding aliases, minimizing the incidence of interlop- 
ers, and providing a reliable estimate of the redshift uncer- 
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Fig. 17. Photometric redshift tests on galaxies in the SWIRE-ESIS field, run with EAZY code (G. Brammer, P. G. van 
Dokkum & P. Coppi, Ap J submitted). Left: all photometric bands; center, all bands excepted / and z\ right: GALEX, / 
and IRAC only. The dashe d lines mark the l t r scat ter in z p h t — z spec . Spectroscopic data belong to Feruglio et al. (sub.), 
Sacchi et al. (in prep.) and lLa Franca et all (|2004h . 



tainties. The EAZY code has the advantage of being very 
user-friendly and specifically designed for multi-wavelength 
surveys. 

Figure [T7] shows the results of these tests, comparing 
photometric and spectroscopic (Fe ruglio et al., sub.; Sacchi 
et al., in prep.; La Franca et al. 2004) redshifts. In the 
left panel, photometric redshift were estimated by using all 
the available UV-optical-IR data (GALEX, ESIS, J, K s , 
IRAC), while the central panel does not account for I and 
z data. Finally the right-hand plot shows photometric red- 
shifts obtained including only GALEX, / band and IRAC 
fluxes (i.e. those data available over the whole ELAIS-S1 
area) . 

When introducing / and z — in combination with all the 
other available data — one would expect a significant im- 
provement in the photometric redshift estimate. Actually, 
these data are effective in constraining photometric red- 
shifts, only above redshift z ~ 0.8, when the Balmer break 
lies in or red ward the R band. Unfortunately, only few spec- 
troscopic redshifts are available in this redshift domain, in 
ELAIS-S1. Nevertheless, one can appreciate some improve- 
ment in the scatter of z p hot — z spec for z > 0.8, where <ta 2 
reduces from 0.36 to 0.25. No significant improvement is 
seen at lower redshifts (ctaz ~ 0.15). 

On the other hand, when considering only GALEX, / 
band and IRAC data, the result reverses. At low redshift 
the match between photometric and spectroscopic z is rea- 
sonably good, and the scatter is not dramatically larger 
than in the previous cases (<taz — 0.17). The code tends to 
systematically underestimate the redshift, but the drift is 
smaller than A(z) ~ 0.05. At redshifts larger than z ~ 0.6, 
the results become rather poor: the scatter increases, and 
the systematic shift is no longer negligible. At z > 0.8 we 
have a systematic Az ~ —0.15. This effect can probably 
be explained with the difficulty to detect galaxies at higher 
redshift with GALEX. An obvious exception is the case of 
IR-peakers, for which the IRAC det ection of 1.6um ste l- 
lar peak constrains th e redshift (e.g. iBerta et al.ir2007alTbl : 
IWeedman et al.ll2006D . 

The fraction of catastrophic outliersj increases from 5% 
to 10%, when using only GALEX, / and IRAC bands. 



8 Catastrophic outliers are denned here as sources having 

Az > 3<7Az. 



Generally speaking, iRowan-Robinson et al.l 
showed that at least 5 different photometric bands are 
needed, in order to obtain a reliable photometric redshift 
estimate, including 2 IRAC fluxes, and (for example) 3 opti- 
cal. ELAIS-S1 will benefit of three additional optical bands 
over the whole field, once the ESIS-WFI data processing 
will be completed and presented in a forthcoming paper. 



7. Summary 

In the era of wide-area cosmological surveys, we have pre- 
sented I and z band wide field imaging of ELAIS-S1, carried 
out with the VIMOS camera on VLT. These observations 
are part of the ESO-Spitzer Imaging extragalactic Survey 
(ESIS) and cover ~ 4 deg 2 in the / band and ~ 1 deg 2 in 
2. The nominal exposure time is 1800s and 3600s in I and 
z respectively, but the coverage is not uniform, due to the 
very complex observing pattern. 

As a result of accurate data processing, nine indepen- 
dent mosaics for the / band and one single, 1 deg 2 wide 
frame in the z band were produced. The final catalogs in- 
clude more than 300000 / band sources and 50000 z band 
objects. Equatorial coordinates of VIMOS objects are con- 
strained within ~ 0.2 r.m.s. arcsec, when compared to the 
GSC 2.2 catalog. Relative astrometric accuracy between the 
nine / band images is as good as ~ 0.1 arcsec r.m.s. 

Photometric uncertainties were estimated through sim- 
ulations, and turn out to be ~ 0.05, 0.20 mag at I = 
21.5, 23.5 and slightly larger in the z band. The survey 
reaches a 90% average completeness at 23.1 and 22.5 (Vega 
magnitudes) in the / and z band respectively. 

The ESIS VIMOS survey provides optical / coun- 
terparts for > 60% (137206/219877) of SWIRE/Spitzer 
sources in the ELAIS-S1 area. 

We have exploited the new I and z band data — in com- 
bination to the other multi-wavelength observations avail- 
able in ELAIS-S1 — to search for high-redshift (z = 1 — 3) 
galaxies. Our selection is based on the detection of the 
1.6/im stellar peak in the IRAC wavelength domain (IR- 
peakers |Berta~eFaI] [2007a). The I band data give a valu- 
able constraint on IR-peakers: requiring (/ — 3.6)ab > 2 
avoids low-redshift [z <J 0.6) interlopers. The availability 
of VIMOS data over ~ 4 deg 2 allows the selection to be 
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extended over the whole SWIRE-ESIS field. A total of 973 
5.8/im-peak galaxies at z — 2 — 3 were found over ~ 4 deg 2 
at the SWIRE 5.8^m 3ct depth (25.8 /xJy). 

The z band d ata were used to perform the BzK 
(|Daddi et al.ll2004T ) selection over ~ 1 deg 2 and verify the 
consistency of the IR-peakers class. Our high-z objects lie 
in the locus of z > 1.4 star forming galaxies. Apparently, 
no passive IR-peakers are detected at SWIRE-ESIS depth, 
over the central ESIS 1 deg 2 field. 

At the SWIRE depth, these galaxies have stellar masses 
M > 10 11 M at z > 2. The identification of a signifi- 
cant number of these rare high-z sources is possible only 
thanks to the wide area probed by SWIRE-ESIS. This sur- 
vey is giving a significant contribution — among others — 
in studying the very massive tail of the stellar mass func- 
tion at z > 2, and constraining the physical properties of its 
constituents (e.g. iBerta et al.ll2007al ; lLonsdale et al.l l2~006) . 

Updates on the ESIS project are available on the web 



page http: //www. astro .unipd. it/esis 
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Appendix A: Treatment of fringing in VIMOS / 
and z band frames 

As mentioned in Section[3l the VIMOS / and z band frames 
suffer a strong fringing that needs to be corrected before 
photometric calibration and mosaicking. Here we describe 
the main causes of fringing and the procedure adopted to 
clean ESIS VIMOS images. 
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A.l. Fringing effects 

At wavelengths shorter than ~7000A, the silicon layer of 
VIMOS CCDs absorbs photons efficiently, close to the sur- 
face of the chip. On the other hand, in the / and z bands 
the coefficient of absorption of silicon falls rapidly with in- 
creasing wavelength. Consequently a significant number of 
low energy photons passes through the silicon, which acts 
as a thin film against the higher index of refraction of the 
substrate: light is reflec ted back and forth bet ween the two 
silicon boundaries (e.g. iMalumuth et al.ll2003l) . 

The crossing and re-crossing waves mutually interfere, 
reinforce or cancel in amplitude, depending on their wave- 
length relative to the film thickness. If the thickness of the 
Si layer varies across the chip, then "fringes" are formed, 
following paths of constant thickness. 

If the CCD is illuminated by mono-chromatic light, with 
a proper wavelength, sharp fringes are produced. If a uni- 
form broad-band source is used, the range of wavelengths 
washes out the appearance of fringes, and their pattern will 
have a reduced amplitude. 

The night sky is dominated by emission lines, particu- 
larly at the red end of the optical domain. The fringe pat- 
tern formed on VIMOS images in the / and z bands is hence 
dominated by the combination of the thin-film interference 
for the different monochromatic sky lines. 

The observed fringe patterns therefore depend not only 
on the silicon thickness across the chip, but also on which 
lines are strong in the night sky. 

Since the composition and conditions of the sky should 
not vary by large amounts across one night, at a first look 
the observed fringe patterns do not vary strongly from 
one image to another. Nevertheless, the amplitude of the 
fringes, and their strength relative to the uniform sky back- 
ground may vary significantly. Any continuum variation in- 
dependent on the line emission does change the uniform 
background relative to the fringe pattern. Some examples 
are excess scattered light from the moon or from the sun 
near twilight. 

In addition to changes in the continuum background, 
also the line strengths may vary, as a consequence of 
changes in the high-altitude atmospheric conditions (e.g. 
temperature, solar- wind particles). In this case, if the line 
ratios change significantly, not only the amplitude of fringes 
changes, but the whole fringe pattern may shift. Even if the 
fringe pattern is mostly stable, second-order changes make 
it very difficult to apply a single pattern to all images taken 
during a full obse rving night (see also the de scription of 
fringing by Elixir. iMagnier fc Cuillandreil2004 ). 



obtained from the science images themselves: by combining 
all frames obtained at different pointings across each ob- 
serving night , it is possible to obtain an improved sky flat 
field fsee lBerta et al.l l2006) containing the fringe informa- 
tion. The procedure adopted to extract the fringe pattern 
is as follows: 

— first of all, mask all astronomical objects from the 
science images, by using the obj masks task in the 
MSCRED IRAF package. 

— very bright stars and extended galaxies produce ex- 
tended halos on the images, which need to be masked 
manually. 

— all science frames are combined together, using object- 
masks. The result is an improved sky flat field, without 
any astronomical source, showing strong fringes. 

— the fringe pattern is then isolated by dividing the im- 
proved sky flat field by a boxcar-median frame. 

The fringe pattern frame is finally subtracted from each 
science image, using an appropriate scaling factor in or- 
der to account for the possible changes in fringes intensity, 
relative to the background, and minimize the residuals. 

If sky conditions were perfectly stable across the whole 
observing night, one could use the same fringe pattern to 
correct all images taken during the night. In practice - 
with the exception of few fortunate nights — sky condi- 
tions are far from being perfectly consistent for the whole 
duration of observations and it is necessary to produce sev- 
eral (three or more) fringe frames per night, splitting it in 
contiguous shorter bits. As a consequence, the number of 
frames used to produce each pattern decreases. 

In the worst cases, it is possible that the frames used to 
produce a given fringe pattern are not pointing in signifi- 
cantly different regions of the sky. This case is particularly 
critical if very bright stars lie in the field of view, because 
there exist regions of the chip where fringes are not sam- 
pled, as a consequence of heavy masking of big haloes. 

An example of the de-fringing procedure is depicted in 
Fig.fJ] where we show one VIMOS frame before (left panel) 
and after (right panel) subtracting the fringe pattern (cen- 
tral panel). 

After fringe removal, it is finally possibl e to apply a 
super-sky-fiat (ss-flat) (see iBerta et al.l [2006T ) in order to 
correct for residual second order inhomogeneities in the 
background. Super-sky-flat frames were produced by com- 
bining all de-fringed science frames and masking astronom- 
ical sources. Again, because of sky variability, several ss-flat 
frames were needed during each night. 



A. 2. De-fringing 



The additive fringing component described in Section IA.1I 
must be removed from all science images, by subtracting a 
template fringe frame, adequately scaled to the amplitude 
of fringes in each image. 

As the light from stars and galaxies is mainly dom- 
inated by continuum emission, with a negligible fraction 
from emission lines, their contribution to fringing is mini- 
mal (see Sect. IA.1|) . and fringes are basically produced by 
the night sky background. 

The ideal fringe pattern frame would be obtained from 
blank sky images taken in the same conditions as science 
frames. The optimal estimation of the fringe pattern is thus 



Berta S., et al.: ESIS-VIMOS 



19 



Field 


OB 


R.A. 


Dec. 


Date 


A.M. 


Field 


OB 


R.A. 


Dec. 


Date 


A.M. 


number 




[hh:mm:ss] 


[»:':"] 


[yyyy-mm-dd] 




number 




[hh:mm:ss] 


[":':"] 


[yyyy-mm-dd] 




Field 1 


A 


00:32:48.00 


-44:24:00.0 


2003-07-08 


1.509 


Field 26 


A 


00:38:26.10 


-43:30:45.0 


2003-11-19 


1.258 




B 


00:32:59.30 


-44:21:36.0 


2003-11-18 


1.063 




B 


00:38:37.24 


-43:28:21.0 


2003-11-18 


1.066 




C 


00:33:10.59 


-44:19:12.0 


2003-07-27 


1.454 




C 


00:38:48.35 


-43:25:57.0 


2003-11-16 


1.090 


Field 2 


A 


00:34:13.80 


-44:24:00.0 


2003-07-08 


1.434 


Field 27 


A 


00:39:50.63 


-43:30:45.0 


2003-11-21 


1.383 




B 


00:34:25.10 


-44:21:36.0 


2003-11-21 


1.239 




B 


00:40:01.76 


-43:28:21.0 


2003-07-10 


1.092 




C 


00:34:36.38 


-44:19:12.0 


2003-07-27 


1.386 




C 


00:40:12.88 


-43:25:57.0 


2003-11-16 


1.165 


Field 3 


A 


00:35:39.59 


-44:24:00.0 


2003-07-08 


1.366 


Field 28 


A 


00:41:15.15 


-43:30:45.0 


2004-08-11 


1.064 




B 


00:35:50.90 


-44:21:36.0 


2003-11-17 


1.270 




B 


00:41:26.29 


-43:28:21.0 


2004-08-21 


1.177 




C 


00:36:02.18 


-44:19:12.0 


2003-07-27 


1.330 




C 


00:41:37.41 


-43:25:57.0 


2004-08-21 


1.068 


Field 4 


A 


00:37:05.39 


-44:24:00.0 


2003-07-08 


1.314 


Field 29 


A 


00:32:48.00 


-43:13:00.0 


2003-08-03 


1.059 




B 


00:37:16.69 


-44:21:36.0 


2003-11-18 


1.085 




B 


00:32:59.08 


-43:10:36.0 


2003-12-19 


1.286 




C 


00:37:27.98 


-44:19:12.0 


2003-07-27 


1.232 




C 


00:33:10.15 


-43:08:12.0 


2004-07-17 


1.201 


Field 5 


A 


00:38:31.19 


-44:24:00.0 


2003-07-08 


1.268 


Field 30 


A 


00:34:12.11 


-43:13:00.0 


2003-12-16 


1.086 




B 


00:38:42.49 


-44:21:36.0 


2003-07-09 


1.268 




B 


00:34:23.19 


-43:10:36.0 


2003-12-21 


1.219 




C 


00:38:53.78 


-44:19:12.0 


2003-11-16 


1.384 




C 


00:34:34.26 


-43:08:12.0 


2004-07-19 


1.104 


Field 6 


A 


00:39:56.99 


-44:24:00.0 


2003-07-08 


1.225 


Field 31 


A 


00:35:36.23 


-43:13:00.0 


2003-12-16 


1.389 




B 


00:40:08.29 


-44:21:36.0 


2003-11-18 


1.065 




B 


00:35:47.31 


-43:10:36.0 


2004-08-21 


1.364 




C 


00:40:19.57 


-44:19:12.0 


2003-07-31 


1.069 




C 


00:35:58.37 


-43:08:12.0 


2004-08-13 


1.389 


Field 7 


A 


00:41:22.78 


-44:24:00.0 


2004-07-12 


1.207 


Field 32 


A 


00:37:00.34 


-43:13:00.0 


2003-12-18 


1.175 




B 


00:41:34.09 


-44:21:36.0 


2004-07-14 


1.074 




B 


00:37:11.42 


-43:10:36.0 


2004-07-22 


1.109 




C 


00:41:45.37 


-44:19:12.0 


2004-08-21 


1.063 




C 


00:37:22.49 


-43:08:12.0 


2004-07-22 


1.067 


Field 8 


A 


00:32:48.00 


-44:06:15.0 


2003-07-08 


1.170 


Field 33 


A 


00:38:24.46 


-43:13:00.0 


2003-07-08 


1.060 




B 


00:32:59.24 


-44:03:51.0 


2003-11-18 


1.064 




B 


00:38:35.54 


-43:10:36.0 


2003-07-10 


1.075 




C 


00:33:10.47 


-44:01:27.0 


2003-07-27 


1.542 




C 


00:38:46.60 


-43:08:12.0 


2003-11-16 


1.242 


Field 9 


A 


00:34:13.37 


-44:06:15.0 


2003-07-08 


1.143 


Field 34 


A 


00:39:48.57 


-43:13:00.0 


2003-11-21 


1.461 




B 


00:34:24.61 


-44:03:51.0 


2003-07-09 


1.144 




B 


00:39:59.65 


-43:10:36.0 


2003-07-10 


1.065 




C 


00:34:35.84 


-44:01:27.0 


2003-07-27 


1.271 




C 


00:40:10.72 


-43:08:12.0 


2003-11-16 


1.133 


Field 10 


A 


00:35:38.73 


-44:06:15.0 


2003-07-08 


1.121 


Field 35 


A 


00:41:12.69 


-43:13:00.0 


2004-08-11 


1.073 




B 


00:35:49.98 


-44:03:51.0 


2003-07-09 


1.116 




B 


00:41:23.77 


-43:10:36.0 


2004-08-20 


1.196 




C 


00:36:01.21 


-44:01:27.0 


2003-07-27 


1.188 




C 


00:41:34.83 


-43:08:12.0 


2004-08-21 


1.144 


Field 11 


A 


00:37:04.10 


-44:06:15.0 


2003-07-08 


1.101 


Field 36 


A 


00:32:48.00 


-42:55:15.0 


2003-08-01 


1.112 




B 


00:37:15.35 


-44:03:51.0 


2003-11-18 


1.099 




B 


00:32:59.03 


-42:52:51.0 


2003-08-03 


1.086 




C 


00:37:26.57 


-44:01:27.0 


2003-07-27 


1.157 




C 


00:33:10.04 


-42:50:27.0 


2003-11-18 


1.117 


Field 12 


A 


00:38:29.47 


-44:06:15.0 


2003-07-08 


1.086 


Field 37 


A 


00:34:11.71 


-42:55:15.0 


2003-11-19 


1.054 




B 


00:38:40.71 


-44:03:51.0 


2003-07-09 


1.085 




B 


00:34:22.74 


-42:52:51.0 


2003-11-18 


1.064 




C 


00:38:51.94 


-44:01:27.0 


2003-11-16 


1.293 




C 


00:34:33.75 


-42:50:27.0 


2003-11-18 


1.139 


Field 13 


A 


00:39:54.83 


-44:06:15.0 


2003-07-08 


1.074 


Field 38 


A 


00:35:35.42 


-42:55:15.0 


2003-11-19 


1.054 




B 


00:40:06.08 


-44:03:51.0 


2003-11-18 


1.060 




B 


00:35:46.45 


-42:52:51.0 


2003-11-19 


1.067 




C 


00:40:17.31 


-44:01:27.0 


2003-11-16 


1.069 




C 


00:35:57.46 


-42:50:27.0 


2003-11-19 


1.058 


Field 14 


A 


00:41:20.20 


-44:06:15.0 


2004-08-18 


1.133 


Field 39 


A 


00:36:59.13 


-42:55:15.0 


2003-11-19 


1.076 




B 


00:41:31.45 


-44:03:51.0 


2004-07-14 


1.065 




B 


00:37:09.70 


-42:52:51.0 


2003-11-16 


1.053 




C 


00:41:42.68 


-44:01:27.0 


2004-08-21 


1.061 




C 


00:37:21.17 


-42:50:27.0 


2003-11-19 


1.081 


Field 15 


A 


00:32:48.00 


-43:48:30.0 


2003-12-17 


1.201 


Field 40 


A 


00:38:22.84 


-42:55:15.0 


2003-11-19 


1.205 




B 


00:32:59.19 


-43:46:06.0 


2003-12-21 


1.138 




B 


00:38:33.86 


-42:52:51.0 


2003-11-16 


1.452 




C 


00:33:10.36 


-43:43:42.0 


2004-07-18 


1.110 




C 


00:38:44.88 


-42:50:27.0 


2003-11-16 


1.534 


Field 16 


A 


00:34:12.94 


-43:48:30.0 


2003-12-16 


1.251 


Field 41 


A 


00:39:46.55 


-42:55:15.0 


2003-11-19 


1.369 




B 


00:34:24.13 


-43:46:06.0 


2004-07-25 


1.058 




B 


00:39:57.57 


-42:52:51.0 


2003-11-18 


1.080 




C 


00:34:35.31 


-43:43:42.0 


2004-07-19 


1.059 




C 


00:40:08.59 


-42:50:27.0 


2003-11-16 


1.104 


Field 17 


A 


00:35:37.89 


-43:48:30.0 


2003-12-17 


1.313 


Field 42 


A 


00:41:10.26 


-42:55:15.0 


2004-07-12 


1.163 




B 


00:35:49.07 


-43:46:06.0 


2004-09-06 


1.064 




B 


00:41:21.28 


-42:52:51.0 


2004-08-20 


1.161 




C 


00:36:00.25 


-43:43:42.0 


2004-09-06 


1.060 




C 


00:41:32.29 


-42:50:27.0 


2004-08-12 


1.148 


Field 18 


A 


00:37:02.83 


-43:48:30.0 


2004-09-06 


1.109 


Field 43 


A 


00:32:48.00 


-44:41:45.0 


2004-06-21 


1.152 




B 


00:37:14.02 


-43:46:06.0 


2004-08-15 


1.164 




B 


00:32:59.36 


-44:39:21.0 


2004-06-23 


1.127 




C 


00:37:25.19 


-43:43:42.0 


2004-08-21 


1.073 




C 


00:33:10.70 


-44:36:57.0 


2004-06-23 


1.106 


Field 19 


A 


00:38:27.77 


-43:48:30.0 


2003-11-19 


1.312 


Field 44 


A 


00:34:14.23 


-44:41:45.0 


2004-06-21 


1.127 




B 


00:38:38.96 


-43:46:06.0 


2003-07-09 


1.063 




B 


00:34:25.59 


-44:39:21.0 


2004-06-26 


1.351 




C 


00:38:50.14 


-43:43:42.0 


2003-11-16 


1.077 




C 


00:34:36.94 


-44:36:57.0 


2004-06-26 


1.290 


Field 20 


A 


00:39:52.72 


-43:48:30.0 


2003-11-21 


1.324 


Field 45 


A 


00:35:40.47 


-44:41:45.0 


2004-06-21 


1.095 




B 


00:40:03.90 


-43:46:06.0 


2003-07-10 


1.140 




B 


00:35:51.83 


-44:39:21.0 


2004-06-26 


1.134 




C 


00:40:15.08 


-43:43:42.0 


2003-11-16 


1.200 




C 


00:36:03.17 


-44:36:57.0 


2004-07-25 


1.082 


Field 21 


A 


00:41:17.66 


-43:48:30.0 


2004-08-20 


1.137 


Field 46 


A 


00:37:06.70 


-44:41:45.0 


2004-06-27 


1.123 




B 


00:41:28.85 


-43:46:06.0 


2004-07-24 


1.069 




B 


00:37:18.06 


-44:39:21.0 


2004-08-15 


1.207 




C 


00:41:40.02 


-43:43:42.0 


2004-08-21 


1.063 




C 


00:37:29.41 


-44:36:57.0 


2004-08-12 


1.098 


Field 22 


A 


00:32:48.00 


-43:30:40.0 


2003-12-17 


1.121 


Field 47 


A 


00:38:32.94 


-44:41:45.0 


2004-06-27 


1.104 




B 


00:32:59.13 


-43:28:21.0 


2003-12-19 


1.457 




B 


00:38:44.30 


-44:39:21.0 


2004-08-12 


1.114 




C 


00:33:10.25 


-43:25:57.0 


2004-07-17 


1.117 




C 


00:38:55.64 


-44:36:57.0 


2004-08-20 


1.097 


Field 23 


A 


00:34:12.53 


-43:30:45.0 


2003-12-16 


1.147 


Field 48 


A 


00:39:59.17 


-44:41:45.0 


2004-06-27 


1.089 




B 


00:34:23.66 


-43:28:21.0 


2004-07-17 


1.061 




B 


00:40:10.53 


-44:39:21.0 


2004-08-20 


1.115 




C 


00:34:34.78 


-43:25:57.0 


2004-07-19 


1.066 




C 


00:40:21.88 


-44:36:57.0 


2004-08-12 


1.133 


Field 24 


A 


00:35:37.05 


-43:30:45.0 


2003-12-18 


1.110 


Field 49 


A 


00:41:25.41 


-44:41:45.0 


2004-08-21 


1.226 




B 


00:35:48.19 


-43:28:21.0 


2004-08-21 


1.142 




B 


00:41:36.77 


-44:39:21.0 


2004-08-21 


1.128 




C 


00:35:59.30 


-43:25:57.0 


2004-08-13 


1.240 




C 


00:41:48.11 


-44:36:57.0 


2004-08-21 


1.069 


Field 25 


A 
B 
C 


00:37:01.58 
00:37:12.71 
00:37:23.83 


-43:30:45.0 
-43:28:21.0 

-43:25:57.0 


2003-12-18 
2004-07-22 
2004-08-13 


1.287 
1.057 
1.066 















Table 4. Log file of / band ESIS VIMOS observations. 
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Field 


OB 


R.A. 


Dec. 


Date 


A.M. 


number 




[hh:mm:ss] 


[":':"] 


[yyyy-mm-dd] 




Field 15 


A 


00:32:48.00 


-43:48:30.0 


2003-12-17 


1.250 




B 


00:32:59.19 


-43:46:06.0 


2003-12-21 


1.174 




C 


00:33:10.36 


-43:43:42.0 


2004-07-18 


1.088 


Field 16 


A 


00:34:12.94 


-43:48:30.0 


2003-12-16 


1.311 




B 


00:34:24.13 


-43:46:06.0 


2004-07-25 


1.063 




C 


00:34:35.31 


-43:43:42.0 


2004-07-19 


1.064 


Field 17 


A 


00:35:37.89 


-43:48:30.0 


2003-12-17 


1.386 




B 


00:35:49.07 


-43:46:06.0 


2004-09-06 


1.076 




C 


00:36:00.25 


-43:43:42.0 


2004-09-06 


1.058 


Field 18 


A 


00:37:02.83 


-43:48:30.0 


2004-09-06 


1.088 




B 


00:37:14.02 


-43:46:06.0 


2004-08-15 


1.130 




C 


00:37:25.29 


-43:43:46.8 


2004-08-21 


1.098 


Field 22 


A 


00:32:48.00 


-43:30:45.0 


2003-12-17 


1.154 




B 


00:32:59.13 


-43:28:21.0 


2003-12-19 


1.561 




C 


00:33:10.25 


-43:25:57.0 


2004-07-17 


1.093 


Field 23 


A 


00:34:12.53 


-43:30:45.0 


2003-12-16 


1.186 




B 


00:34:23.66 


-43:28:21.0 


2004-07-17 


1.057 




C 


00:34:34.78 


-43:25:57.0 


2004-07-19 


1.058 


Field 24 


A 


00:35:37.05 


-43:30:45.0 


2003-12-18 


1.138 




B 


00:35:48.19 


-43:28:21.0 


2004-08-21 


1.181 




C 


00:35:59.30 


-43:25:57.0 


2004-08-13 


1.191 


Field 25 


A 


00:37:01.58 


-43:30:45.0 


2003-12-18 


1.372 




B 


00:37:12.71 


-43:28:21.0 


2004-07-22 


1.061 




C 


00:37:23.83 


-43:25:57.0 


2004-08-13 


1.059 


Field 29 


A 


00:32:48.00 


-43:13:00.0 


2003-08-03 


1.055 




B 


00:32:59.08 


-43:10:36.0 


2003-12-19 


1.355 




C 


00:33:10.15 


-43:08:12.0 


2004-07-17 


1.159 


Field 30 


A 


00:34:12.11 


-43:13:00.0 


2003-12-16 


1.107 




B 


00:34:23.19 


-43:10:36.0 


2003-12-21 


1.273 




C 


00:34:34.26 


-43:08:12.0 


2004-07-19 


1.082 


Field 31 


A 


00:35:36.23 


-43:13:00.0 


2003-12-16 


1.478 




B 


00:35:47.31 


-43:10:36.0 


2004-08-21 


1.293 




C 


00:35:58.37 


-43:08:12.0 


2004-08-13 


1.314 


Field 32 


A 


00:37:00.34 


-43:13:00.0 


2003-12-18 


1.219 




B 


00:37:11.42 


-43:10:36.0 


2004-07-22 


1.086 




C 


00:37:22.49 


-43:08:12.0 


2004-07-22 


1.058 



Table 5. Log file of z band ESIS VIMOS observations. 



